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Single photons are a vital resource for photonic quantum information processing. However, even 
state-of-the-art single photon sources based on photon-pair generation and heralding detection have 
only a low probability of delivering a single photon when one is requested. We analyse a scheme that 
uses a switched fibre delay loop to increase the delivery probability per time bin of single photons 
from heralded sources. We show that, for realistic experimental parameters, combining the output 
of up to 15 pulses can yield a performance improvement of a factor of 10. We consider the future 
performance of this scheme with likely component improvements. 


Heralded single photon sources based on photon-pair 
generation in nonlinear media are at the forefront of de¬ 
velopments in photonic quantum technologies [l[. Pho¬ 
tons are created in pairs, usually by a pulsed pump laser, 
allowing the delivery of a single photon in a well-defined 
time bin to be conditioned on the detection of its twin Q . 
The simplicity of this approach is attractive, but, because 
pair generation is spontaneous, the generation probabil¬ 
ity must be kept low in order to limit contamination 
from higher-order photon-pair components [^. Com¬ 
bined with the effects of loss, the probability of delivering 
a heralded photon from any individual pump laser pulse 
is typically less than 1%. 

Realistic quantum information protocols require man’ 


if 


single photon sources to be operated simultaneously 
The probability of this occurring, P, for N sources be¬ 
comes vanishingly small if the individual probability of 
each source producing a photon p <^\ because P = p^. 
One method of increasing p is by multiplexing - actively 
combining the output from several generation modes us¬ 
ing delay, feed-forward, and a fast switch Com¬ 

bining the output of pair generation in several separate 
sources - spatial multiplexing - brings significant perfor¬ 
mance benefits, but the resource overhead in nonlinear 
media and switches is high ii- 

Alternatively, one can combine the output from dif¬ 
ferent temporal modes of the same source, bringing an 
overhead cost only in time and not in additional physical 
resources In this letter, we present a novel mul¬ 

tiplexing scheme in the temporal domain based on a fibre 
loop that uses only a single optical switch and delay line. 
We numerically evaluate the performance of the scheme 
in the presence of imperfect detection, switch loss and 
attenuation in the delay line by extending the methods 
we presented in [l^ and consider the realistic limits to 
source performance. 

Fig-IDshows the proposed setup. A non-linear medium 
(^(2) crystal or optical fibre) is pumped by a pulsed 
laser. The generated signal-idler pair is first split in wave¬ 
length, the signal photon is sent to a detector to herald 
the presence of the remaining idler photon. In this work, 
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FIG. 1. Schematic of temporal multiplexing scheme. See text 
for details. 


we consider a detector that has photon number resolving 
(PNR) capability, though the analysis could be equally 
well applied to a binary detector. The heralded idler pho¬ 
ton is coupled into fibre and subjected to a fixed delay to 
give sufficient time for detection and feed-forward. Fol¬ 
lowing the fixed delay a 2 x 2 optical switch is used to 
route the photon either to the output or into a fibre stor¬ 
age loop with a propagation delay of exactly one period 
of the the laser pulse train. Hence the probability am¬ 
plitude of photons held in the storage loop will overlap 
with those of heralded idler photons from any subsequent 
laser pulse. 

The multiplexing scheme is described as follows. A 
train of m laser pulses, labelled from t = 1 to t = m in 
order of arrival, enters the nonlinear photon-pair gener¬ 
ation medium. If a pair is generated on the first pulse 
and the heralding detector fires, the switch is set to the 
crossed state so that the corresponding idler photon is 
routed into the the storage loop. On the next pulse if the 
heralding detector fires again the switch is once again set 
to the crossed state, simultaneously rejecting the photon 
from the previous pulse and storing the new photon in 
the loop. If there is no heralding event, then the switch 
remains closed and photon in the storage loop completes 
another pass through the switch. This process is repeated 
up to and including the m’th pulse and we refer to m as 
the multiplexing depth. 

On the m’th pulse, if there is no successful heralding 
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event then the switch routes the stored photon amplitude 
into the optical output. However, if there is a successful 
heralding event on the final pulse the switch allows the 
newly-generated photon straight through to the output. 
Hence the photon leaving the output following the m’th 
pulse may have passed through the switch and storage 
loop anywhere from t = 1 to t = m times and therefore 
accrued a vastly different amount of loss. At low values 
of t a photon must make more passes through the loop 
to be used after the m’th pulse and so its contribution to 
the overall probability of successfully delivering a single 
photon from the scheme is lower compared to photons 
generated on later pulses where t —>■ m. Furthermore, the 
loss of the switch and storage loop have a large effect on 
the probability of successfully delivering a single photon. 

To evaluate the performance of the proposed system 
we follow a similar approach as outlined in . We first 
assume that the photon-pair source has been engineered 
to generate signal and idler photons in only two spatio- 
temporal modes such that the modes are only correlated 
in photon number Q. This allows us to describe the 
probability amplitudes of the state-vector, 

|d^)=Q'o|0s,0i)-|-ai|ls,li)+fl2|2s,2i)-|---' , (1) 

with thermal statistics: 

= ( 2 ) 

[n + l) \n + l/ 

where n is the mean photon number per pulse. 

A simple result can be obtained by only considering 
the first non-zero term in the state-vector expansion in 
Eq. [T] The overall probability of successfully delivering 
a heralded single photon on the m’th pulse, p(success), 
is given by 

m 

p(success) = 1 - n “ Pth{l)VdVL) : (3) 

t=i 


which are dependent on the number of photons rig de¬ 
tected by the heralding detector. By working with a 
PNR detector all detection results where rig ^ 1 can be 
ignored. To extract the probability of successfully deliv¬ 
ering a single photon from the t’th pulse we multiplied 
Pii^s = l,t) by the probability of achieving a heralding 
detection result of rig = 1 and calculated the overlap with 
a pure single photon Fock state: 


^(success, t) = p{ns 


Tr{pi{ns = l,t)} ' 


(4) 


The overall probability of successfully delivering a single 
photon from a pulse train divided into different tempo¬ 
ral bins, each containing m pulses is found through the 
product 


m 

^(success, m) = 1 — (1 — ^(success, t)). (5) 

t=i 

The contribution to noise, defined as a successful herald¬ 
ing detection followed by the delivery of an idler state 
containing more than one photon, was found through. 


p(noise, m) 



OO 


E 


{n\pi{ns = l,t)\n) 
Tr{pi(ng = l,t)} 


When Pi is properly normalized such that Tr{pi} 
this reduces to: 


( 6 ) 


= 1, 


m / 1 

p(noise,m) = 1 “ 

i—1 \n—0 



(7) 


it is nevertheless important to note that the accuracy 
of both ^(success, m) and p(noise, m) will depend on the 
value of n at which the calculation is truncated due to 
the effect this has on the normalization of pi. We define 
a corresponding signal-to-noise ratio. 


where ijd and tjl are the detector efficiency and lumped 
efficiency of one pass of the switch and storage loop. 
However, this expression is only valid at low values of 
h. In order to describe the behaviour of the source 
more faithfully at higher mean photon numbers we have 
analysed the effects of higher photon-number terms from 
the state vector. To do so we applied the formalism 
developed in 1^ in which we determined the effect of 
concatenated loss in spatially multiplexed single photon 
sources. We used this to calculate the reduced density 
matrix, pi (rig , t ), describing the heralded idler state given 
a heralding detection of rig photons on the t’th pulse and 
delayed until the m’th pulse. We modelled the effects of 
loss by applying a beam splitter transform in which the 
transmission coefficient is given by the lumped efficiency 
of t passes through the switch and delay line. 

For each pulse t there exists a set of normalised density 
matrices {;di(ng,<) : t,ng G N} with Tr{;di(ng,<)} = 1, 


SNR = 

p(noise, m) 


( 8 ) 


that allows us to make direct comparisons of systems 
with disparate characteristics. Using these techniques, 
we have simulated numerically the effects of higher-order 
photon-pair components on a temporally-multiplexed 
source to enable the optimization of real-world systems 
in the presence of loss. 

Figure [2] shows the performance of the proposed mul¬ 
tiplexing scheme with rja = 0.7 and pL = 0.8 for differ¬ 
ent multiplexing depths m. As the multiplexing depth 
is increased the resulting probability of successfully de¬ 
livering a single photon from the output increases. The 
corresponding SNR is shown in Fig. |2(b)[ To maintain 
a high SNR the mean photon number is constrained to 
be low in order to reduce the contribution from multi¬ 
photon components in the heralded density matrix. By 
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FIG. 2. a.) Variation in p(success) with increasing n for mul¬ 
tiplexing depths of m = 1 (bine), 3 (green), 5 (red) and 15 
(purple) pulses, b.) Corresponding signal-to-noise ratio, c.) 
Difference (Ap) between p(success) calculated using Eq|3]and 
the full numerical calculation with up to 5 photon number 
components. Need to label the lines somehow d.) The in¬ 
crease in p(success) with increasing multiplexing depth at a 
fixed n (green squares) and fixed SNR = 100 (blue circles). 
Triangles show the comparison with a spatial multiplexing 
scheme for fixed h and fixed SNR in teal and red respectively. 
Simulations carried out with a detector efficiency r]d = 0.7 
and lumped loop efficiency r]L = 0.8. 


multiplexing, the overall probability of successfully deliv¬ 
ering a single photon at the output is increased at fixed 
SNR, as seen in Fig. |2(d)[ 

It is insightful to find the limits of validity to the sim¬ 
ple analytic result of eq. [3] and thence to set the point 
at which our numerical calculations can safely be trun¬ 
cated. The difference between the analytic result and the 
numerical density matrix method truncated at different 
numbers of photon pairs is plotted in Fig. 2(c) We can 
see that for larger mean photon numbers it is essential 
to include higher-order terms in the density matrix, how¬ 
ever, the benefit in using any amplitudes above 5 photon 
pairs quickly diminishes. As a result, all subsequent cal¬ 
culations were made by truncating the state vector at 
n = 5. 


Figure [2(d)] shows the increase in p(success) with mul¬ 
tiplexing depth m. Results from our previous simulations 
of a spatial multiplexing scheme are shown for compari¬ 
son |15l | . Our numerical approach allows a direct compar¬ 
ison between multiplexed devices at a fixed SNR of 100. 


It can be seen that the temporal multiplexing loop yields 
an increase in p(success) commensurate with spatial mul¬ 
tiplexing schemes of approximately the same depth, but 
with a large reduction in the resource overhead required 
for implementation. In the case of temporal multiplex¬ 
ing the exponentially diminishing gain in p(success) with 
increasing passes through the storage loop leads to the 
improvement saturating at larger values of h as shown in 
Fig. |2(d)] For small values of h < 0.01, as the multiplex¬ 
ing depth increases, p(success) tends to a constant value 
given by 


lim {^(success, m)} = 

m—)-oo 




I-tIlJ 


) ^(success, 1). (9) 


This is in contrast to spatially-multiplexed sources, in 
which the probability of success continues to increase 
with multiplexing depth over the range studied here; 
nevertheless, as shown in [l5| . even in the spatial case 
p(success) will begin to decrease at a value dependent on 
the concatenated switch loss. 

At a fixed SNR of 100, a storage loop and switch with 
total loss rjL = 0.8 enables a factor of 6 increase in per¬ 
formance at a multiplexing depth of m = 8 pulses. This 
requires only one photon-pair source, one detector, and 
one switch. To achieve a similar increase in performance 
from a spatially-multiplexed scheme would require vastly 
more resources: 8 sources, 8 detectors, and 7 switches and 
delay lines. Not only is building such a network compli¬ 
cated and expensive, but the sources and delays must 
be carefully matched to ensure that the photons output 
are in identical pure states. The resource scaling for the 
schemes discussed is outlined in Table ID 

We note that after temporal multiplexing has taken 
place the number of time bins in which a photon can be 
delivered has been reduced by a factor of m. Therefore, 
in terms of number of heralded single photons delivered in 
one mode per second, this scheme may be outperformed 
by non-switched photon-pair sources pumped by a high- 
repetition-rate laser [l6|. Nevertheless, if our goal is to 
deliver several single photons from independent sources 
simultaneously the scheme will yield huge improvements 
over current source architectures. This is summarized 
in terms of the waiting time to deliver N single photons 
from N independent sources shown in Fig. [3] 

It is clear that to construct the highest performance 
single photon sources the loss of all the components be¬ 
tween the point of generation and the output must be 
minimised. Switch loss is critical in this endeavour. 
Notwithstanding this, our analysis shows that, in the 
light of the inevitable imperfections in currently-available 
components, meaningful gains in source performance can 
be readily achieved by implementing this scheme. 

Our method of analysis at a fixed SNR enables the 
effects of detector inefficiency and switch loss to be 
straightforwardly compared. Fig. 0] shows the behaviour 
of p(success) with varying detector and switch efficiencies 
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TABLE I. Multiplexing scheme performance comparison. Improvement calculated for multiplexing depth m = 8, relative to 
single source with a heralding detector of efficiency r/d = 0.7 and switch efficiency rjL = 0.8. 


Scheme 

Sources 

Heralding Detectors 

Switches 

Rep. Rate 

SNR 

n 

Improvement 

Temporal 

1 

1 

1 

Rp/d 

100 

- 

6.06 

Spatial 

2d 

2d 

2^* - 1 

Rp 

100 

- 

7.40 

Temporal 

1 

1 

1 

Rpjd 

- 

0.01 

3.32 

Spatial 

2d 

2d 

2^* - 1 

Rp 

- 

0.01 

4.03 



No. of Independent Photons 


FIG. 3. The waiting time to deliver N single photons from 
N independent sources at fixed SNR = 100. Blue (solid) - 
Shows the waiting time for a standard single-photon source 
pumped at 80MHz. Green (circles) - A multiplexing scheme 
with detector efficiency rjd ~ 0.7, switch efficiency r/L = 0.8 
and a multiplexing depth of 15. Red (squares) - A future 
multiplexing scheme with detector efficiency rjd = 0.98, switch 
efficiency rjL ~ 0.95 and a multiplexing depth of 15. For large 
numbers of independent photons, these temporal multiplexing 
schemes exhibit a distinct advantage over traditional single 
photon sources. 


at a fixed SNR = 100 at a fixed multiplexing depth of 
10 . 


We see in Fig. 4(a) that at low values of rjd the de¬ 


tector has limited ability to discriminate between pulses 
on which a single pair or multi-pair generation event has 
occurred. The source must then be operated at very low 
mean photon number in order to maintain the signal to 
noise, severely limiting the success probability. As the 
detector efficiency increases, the probability of success 
grows rapidly due to the larger value of mean photon 
number that can be accessed; as 77 ^ —>■ 1 the source can 
be operated where the one-pair component of the thermal 
distribution of photon number peaks. In this regime we 
are only limited by loss in the switch and the attenuation 





FIG. 4. a.) Variation ofp(success) with detector efficiency, rjd, 
for a fixed switch efficiency {rjL ~ 0.8) at fixed SNR = 100. 
b.) Variation of p(success) with switch efficiency, r/s, for a 
fixed detector efficiency {rid = 0.7) at fixed SNR = 100. All 
points correspond to multiplexing depth of 10. c.) Realistic 
future device with rjd = 0.98 and rjL = 0.95, multiplexing 
depth 1 (solid), 5 (dashed), 10 (dotted) and 15 (dot-dash) 
which approaches deterministic operation, d.) Probability of 
success surface plot showing the overall effect of component 
efficiency for a multiplexing depth of 10 at a fixed SNR = 100. 


of the fibre in the storage loop. In the case of Fig. |4(b)[ 
provided the switch has an efficiency greater than 0.5 
then multiplexing is beneficial, however if the switch ef¬ 
ficiency is not much greater than this limit, multiplexing 
over a large number of pulses will be not be worthwhile 
as the switch loss will throttle any improvement. Con¬ 
versely, if the efficiency of the switch is high then the 
scheme will be limited by the efficiency of the heralding 
detector; the source should be pumped at a mean photon 
number that is low enough to limit the contribution of 
multi-photon events to the output, but multiplexed over 
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a large number of pulses to raise the overall probability of 
success. Fig. 4(c) shows the potential for a future realis¬ 
tic device, with optimised detector and switch efficiencies 
to deliver a near deterministic single photon source using 
this novel architecture. 

In conclusion, we have shown that a temporal loop 
multiplexing scheme is an attractive prospect for increas¬ 
ing the performance of heralded single photon sources 
whilst minimising the number of components and the 
complexity of the complete device. Our analysis shows 
that to realise the benefits of multiplexing, the pump 
must be correctly tailored to the loss of the device. 
Therefore, numerical simulations such as those carried 
out above are key to extracting the maximum perfor¬ 
mance of these systems. By working in an all-fibre archi¬ 
tecture loss can be minimised while ensuring straightfor¬ 
ward integration. 
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